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ABSTRACT 


Longitudinal and transverse horizontal wind velocity 
variability is parameterized based on data collected for the 
Bureau of Land Management (BIM) between 1980 and 1982. The 
BLM data sets were used to derive horizontal wind velocity 
components and standard deviations of these components in 
the longitudinal (downwind; u,0,,) and transverse (crosswind; 
v,%,) directions for various time averages. The standard 
deviations were plotted versus wind speed so that trends in 
the change of values could be determined. A methodology was 
used to develop a theoretical algorithm to account for 
turbulence production mechanisms due to shear, buoyancy and 
mesoscale processes. Empirical fits were determined for 
both the u and v wind directions over four different time 
averages. The results were related to puff growth in the x, 
y directions and form the basis for improving a puff 


dispersion model for near-coast overwater regions. 
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I. INTRODUCTION 


Modeling atmospheric transport within the marine 
atmospheric boundary layer (MABL) has progressed in recent 
years. The uses of such models are numerous, but for the 
purposes of this thesis, they were oriented towards the 
diffusion of airborne contaminants for use in industrial 
applications such as gaseous releases from offshore oil 
Grilling platforms, and Department of Defense uses such as 
chemical-warfare defense. To develop ach models, people 
have taken into account certain geophysical, thermodynamic, 
dynamic processes and developed mathematical algorithms to 
describe the transport of airborne particles. Previous 
works by Hojstrup (1982) and Panofsky and Dutton (1984) 
provided expressions for atmospheric turbulence intensity as 
a function of shear and buoyancy production mechanisms. 
Schacher, et al. (1982) performed a series of experiments 
for the Minerals Management Service (formerly the Bureau of 
Land Management (BLM)) to obtain empirical data to evaluate 
the existing models and to parameterize Environmental 
Protection Agency (EPA) approved Gaussian and trajectory 
dispersion models specifically for the Southern California 
coastal regions. Schacher et al. (1986) further utilized 
these data to parameterize horizontal wind direction 


variability with three fitting parameters to represent 


shear, buoyancy and mesoscale production mechanisms. The 
purpose of those efforts was to improve the prediction of 
gaseous dispersion within the MABL. 

The purpose of this thesis is to parameterize 
longitudinal and transverse horizontal wind velocity 
variability on the basis of the BLM experiments and relate 
them to longitudinal and transverse puff release growth. The 
BLM data sets were used to derive horizontal wind velocity 
components and standard deviations of these components in 
the longitudinal (downwind; u,0,,) and transverse (crosswind; 
V,0yv) directions for various time averages. The standard 
deviations were plotted versus wind speed so that trends in 
the change of values could be determined. A methodology 
partially based on Hojstrup (1982), and Schacher et al. 
(1986) was used to develop a theoretical algorithm which 
accounted for turbulence production mechanisms due to shear, 
buoyancy and mesoscale processes. Empirical fits were 
determined for both the u and v wind directions over four 
different time averages. The results were related to puff 
growth in the x, y directions and form the basis for 
improving a puff dispersion model for near-coast overwater 


regions. 


II. DATA ACOUISITION 


A. GENERAL DESCRIPTIONS 
Data used in this work were collected during the 
California Coastal Offshore Transport and Diffusion 
Experiments in 1980, 1981, and 1982. Four experiments were 
performed, two near Ventura, CA in the Santa Barbara 
Channel, and two near Pismo Beach, CA, in an open coastal 
area. These sites were chosen to be representative of areas 
encountered along the coast. The Environmental Physics 
Group of the Naval Postgraduate School (NPS) performed 
measurements aboard the research ship R/V Acania and was 
responsible for collecting the overwater meteorological data 
for all four experiments. Studies were performed such that 
both the winter and summer seasons could be examined. 
Complete descriptions of the two phases of the experiment 
were presented by Schacher, et al. (1982). To conform with 
previous reports, the experiments of Phase I are referred to 
as BLM 1 and 2, and the experiments of Phase II as BLM 3 and 
4. To make this thesis self-contained, portions of the 
Schacher, et al. (1982) report are included in this thesis. 

1. Description of Experiments 

The experimental approach was as follows: SF6 gas 
Was released within the outer continental shelf area, 


outside of three nautical miles (n mi) from shore. During 


Phase I near Ventura the ship was approximately five n mi 
from shore, whereas during Phase II off Pismo Beach the 
distance was three n mi. The plume location was determined 
by continuous analyzers in aircraft and ground vehicles, by 
grab samples from a boat and on land, and by stationary 
one-hour average samplers on land. The release ship was 
equipped with a complete set of meteorological instrumenta- 
tion, including radiosondes, to determine overwater 
conditions. An aircraft was used to get soundings of mean 
meteorological parameter profiles. Onshore instrumentation 
included fixed and tetroon borne sensors and a Doppler 
acoustic sounder for determining wind profiles. Charts 
showing the locations of the ship, aircraft trajectories, 
and ground level sampling during Phase II are in Fig. 2-1. 
The sequence of events on a sampling day was 
approximately as follows: Throughout the early morning the 
ship reported winds to the shore command station. These 
data and shore wind informations were used for initial 
positioning of the ship. Continuous monitoring of the wind 
was done to determine when the sea breeze had become well 
established and to position the ship so that at the 
shoreline the plume would intersect the center of the fixed 
sampler array. The ship was anchored at a fixed position 
since movement of the ship during a release would introduce 
apparent meander to the plume, contaminating the test 


results. The tracer gas release usually began between 1100 
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and 1300 hours. Mobile sampling began about one-half hour 
after the start of the release. A complete experiment 
lasted 6-8 hours. The tracer gas was monitored continuously 
at the source to maintain a constant flow rate. 

For all experiments, an aircraft carrying a 
continuous SF6 analyzer made near shoreline transects of the 
plume at several elevations. The instrument provided 
readings of the instantaneous concentrations of SF6 in the 
ambient air as a function of position and time. This 
allowed the plume dimensions to be defined in both the 
horizontal and vertical directions. Ground level transects 
were made by a similarly equipped van, operating at the 
shoreline and/or inland. Ground level plume concentrations 
were determined by placing fixed one-hour average collectors 
along one or more fixed arrays parallel to the shoreline. 
These samplers were placed close enough together af that 
several would be within the narrowest plume expected. The 
array was wide enough so that the plume would be within its 
extent even if considerable meander occurred. This 


necessitated the use of a large number of samplers. 


B. ENVIRONMENTAL CONDITIONS 
1. Description of Localities 
In this section the general nature of the Pismo 
Beach locality and climatology is described to identify the 
expected meteorological conditions. Descriptions of the 


synoptic and local conditions pertaining to the Phase II 


experiments are included. Locality descriptions and local 
conditions during the Phase I ee nants are omitted since 
the data sets from Phase I were not utilized in this work. 
A chart of the geographical area surrounding Pismo Beach is 
shown in Fig. 2-2. 

The test locations were chosen (1) because they are 
representative of important types of coastal areas, and (2) 
because both are candidates for new or increased outer 
continental shelf oil development. Thus, the areas afford 
the opportunity to investigate transport under differing 
meteorological conditions and also satisfy the needs for BLM 
modeling for regulatory purposes. 

a. Pismo Beach, CA (Phase IT) 

Pismo Beach is approximately 50 miles north of 
Pt. Conception, ina fairly open coastal area. Pt. Buchon, 
immediately to the north, has 1000 to 2000 foot high hills 
projecting some five miles out to sea. The point influences 
the local flow somewhat but the influence appears to be 
slight. The immediate inland hills are low giving a weaker 
land-sea breeze cycle than near Ventura. The experiments 
were carried out at the mouth of the Santa Maria valley, 
which steers the local flow slightly. The entrance to the 
valley at the beach is approximately eight miles wide and 
the immediate hills on each side of the valley are only one 
to two hundred feet high, so their effect is small. The 


area is representative of an open California coastal region 
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Fig. 2-2 Pismo Beach, California Area 


where air masS movement is controlled by the synoptic 
pressure gradient, giving predominantly northwest flow with 


a long over-water fetch, and by the land-sea breeze cycle. 


2. Synoptic-Scale Seasonal Climatology of the 
California Coast 


The synoptic climatology is the same for the two 
test areas. General descriptions, according to seasons, for 
both areas are as follows. 

a. Summer 

The North Pacific semipermanent subtropical high 
lies to the west of the area and controls the synoptic-scale 
Paso . Clockwise flow around the high produces 
northwesterlies along much of the coast, with the local 
sea-breeze turning the wind more westerly. The general 
onshore flow is aided by the inland thermal trough which is 
created by overland heating. Strong subsidence creates the 
prevalent capping inversion and the occasional passage of 
weak upper-level troughs will dissipate or lift the 
inversion for periods of 12-24 hours. 

be Fadl 

The building of high pressure in the Great Basin 
causes frequent Santa Ana conditions. Santa Ana conditions 
are characterized by high pressure areas located over the 
southern California deserts and lower pressure situated over 
the southern California coastal regime. The pressure 
differential between the desert and coast causes very 


strong, dary, low-level easterly winds to prevail. The 


intensity of these winds is exacerbated by being funneled 
through the numerous canyons, in effect resulting in a hot, 
arr jet. This jet produces much vertical, shear induced 
turbulence in the coastal regime MABL. The pattern of 
storms and upper-level westerlies moves further south 
breaking up See summer pattern. Frontal passage becomes 
more frequent than in the summer and the subtropical high 
becomes displaced or shrinks, resulting in a break up of the 
coastal marine inversion. 

c. Winter 

Frontal passage becomes much more frequent and 

strong surface westerlies often follow the passage. Santa 
Ana winds can still occur when the sea-level pressure in the 
Great Basin becomes sufficiently high. Also, the Pacific 
High and capping inversion can reform between frontal 
passage pecureenceas 

a. Spring 

As the storm pattern moves north, the Pacific 
High again becomes the dominant feature. Cold lows pass 
frequently, followed by strong westerlies. 
3. Wind Climatology 

Wind climatologies are useful in determining 
expected conditions and for assessing whether observed 
conditions are typical. It is not possible to use the 
climatology to accurately predict local conditions on a day- 


by-day basis but seasonal patterns are quite reproduceable. 
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In coastal areas, conditions differ from location to 
location so that site specific climatologies are needed. 
Climatological data are presented for Vandenberg Air Force 
Base (AFB), CA (Det 30, 2WS, 1982) which is near Pismo 
Beach. Vandenberg AFB is approximately 20 miles from the 
Pismo Beach experiment area. No closer coastal climatology 
is available. The data set was obtained at meteorological 
stations which are two to three miles inland. 

Monthly wind averages for Pt. Mugu, CA, including 
the number of days of occurrence of Santa Ana conditions, 
are presented in Table 2-1. Santa Ana conditions are 
widespread so that these data would be approximately correct 
for Vandenberg AFB also. 

Surface wind roses for three-month periods for 
Vandenberg are shown in Fig. 2-3. Wind speed is indicated 
by the width of each "vector," wind direction by the angle, 
and frequency of occurrence by the length. The numbers in 
Poewmmincenose circle are the percentage of the time the 
wind is < 3 knots. The wind speed averages for Vandenberg 
AFB are always less than 10 knots. 

4. Test Period Weather Description 

The following is a description of the synoptic and 
local conditions during the Phase II test periods. The 
synoptic conditions were derived from daily weather maps 
(weekly series), published by the National Oceanic and 


Atmospheric Administration. The description of local 
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TABLE 2-1 


MONTHLY AVERAGES OF THE MOST FREQUENTLY OBSERVED WIND 
DIRECTION, PERCENTAGE OF TIME THE WIND SPEED IS GREATER 
THAN 21 KNOTS, AND THE NUMBER OF DAYS OF SANTA ANA 

THE MAXIMUM NUMBER IS THE MAXIMUM 
OBSERVED OVER A TEN-YEAR PERIOD. 


WINDS PER MONTH. 


Most 
Frequent 

Wind Speed 
Month Direc. (Kt) _ 
JAN NE 10 
FEB W 9 
MAR W 10 
APR W 10 
MAY W 9 
JUN W 8 
JUL W 7 
AUG W 8 
SEP W a 
OCT W 7 
NOV NE 7 
DEC N 5 
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Fig. 2-3 Climatological Wind Roses for Vandenberg AFB 
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conditions is based on the overwater meteorological data 
obtained on the R/V Acania which was operated by the 
Department of Oceanography of the Naval Postgraduate School. 
a. Synoptic Descriptions 
(1) BIM-3 (December 1981). Synoptic-scale 
features and associated west coast flow patterns were 
typical for this time of the year. An upper-air north-south 
ridgeline over the western states was the dominant feature 
and led to generally weak surface pressure gradients off the 
southern California coast. The Mexican thermal low and 
afternoon sea breeze determined the flow associated with the 
ridge's presence. Also typical for the time of year was the 
passage of a fast moving upper wave, and associated 
precipitation and moderate northwest winds. Another wave 
waS approaching at the end of the period. More detailed 
descriptions of the synoptic-scale features and resulting 
coastal flow pattern follows. 
On 7 December a 500 mb ridgeline extended 
North-South from eastern British Columbia to southern 
California. It had a slow eastward progression in advance 
of an approaching upper-level trough extending southward 
from a closed low centered over the Gulf of Alaska. Coastal 
winds on 8 December were easterly during most of the day due 
to the unusual location of the Mexican thermal trough. 
Light westerly winds occurred from 1000 to 1600 in 


conjunction with the local sea-breeZe. 
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The approaching upper-level trough crossed 
the west coast on 10 December and a surface front passed the 
experimental area late on the same day. An extensive 
precipitation area existed along the west coast from 
Southern California to Washington state. Coastal winds 
progressed from southerly on 9 December to northerly on 10 
December with the frontal passage. The northerly gradient 
flow behind the front combined with the afternoon sea breeze 
led to a maximum onshore wind of 15 kts during the afternoon 
of 11 December. 

A weak North-South ridgeline was 
re-established over the West Coast on 11 December and 
existed through 14 December. A fast moving upper level 
short wave (trough) moved through the weak ridge and crossed 
the West Coast in the vicinity of northern Washington on 15 
December. The associated surface front reached northern 
California but did not affect the experimental area. During 
the 12 to 16 December period the coastal wind directions and 
Speeds exhibited flow associated with a Mexican thermal 
trough and the afternoon sea breeze. The winds were east to 
northeast except for the afternoon (1100-2400) when onshore 
flow (northwest) occurred with speeds of 10-15 kt. 

The upper-level ridge intensified on 16 
December and the associated large surface high region 
extended from eastern British Columbia to Nevada. The 


increased gradient on the western side of the surface high 
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led to general easterly winds on 17 December with no 
discernable influence from the sea-breeze effects. 

(2) BILM-4 (June 1982). The synoptic-scale 
conditions and resulting precipitation and coastal wind 
regimes were atypical for the early summer season. The 
Mexican thermal trough should dominate this region, with 
resulting light coastal winds influenced by the sea-breeze 
during this period. Two upper-level troughs passed over the 
west coast during the period. The first (22 June) was a 
fast moving short wave and the second (28-30 June) was a 
deep system associated with a closed low at 500 mb, which 
became nearly stationary over central California. Both 
systems had considerable north-south extent which led to the 
southern California surface pressure patterns reflecting 
their passage. This resulted in a greater than normal 
offshore pressure gradient and a fairly steady onshore wind, 
lacking the usual strong land-sea breeze cycle. Hence, 
strong onshore winds occurred. 

During the 21-23 June period, an upper- 
level trough was moving from off the west coast into the 
mountain states. A surface trough extended from western 
British Columbia into northern California which was an 
intensification of the northern extension of the Mexican 
thermal trough. 

During the 24-26 June period, a more 


intense upper-level trough was developing off the west 
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coast. A surface trough line extended from upper Mexico to 
Washington and the offshore pressure gradient was moderate. 
By 25 June, a closed 500 mb low had formed west of the 
Oregon-Washington coastline. The trough and low intensified 
as they progressed slowly eastward. Precipitation was 
observed along the northern California coast on both 25 and 
26 June and a cold front existed over northern California on 
26 June. 

During the 27 June--1 July period, the 
upper-level trough and closed low moved across the west 
coast. The closed low moved south-eastward and was centered 
over central California on 30 June. Widespread 
precipitation occurred over the central and northern 
California coastal regions on 27 and 28 June and extended 
southward into the experimental area during the 29-30 June 
period. A closed surface low and associated frontal systems 
formed over northern Utah on 28 June. Baeause of the 
offshore pressure gradient, winds remained southwest to 
northwest during this period. They were maximum on 27 June 
(15-20 kt) and decreased gradually, to a 5-10 kt range, on 
the remaining days as the above systems moved across the 


coast and inland. 


C. EQUIPMENT 


1. R/V Acania 


The R/V Acania served as the platform for the 


offshore release of the SF6 tracer gas, and provided 
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continuous measurement of several critical meteorological 
parameters. These measurements were performed to document 
atmospheric transport and stability conditions of the 
overwater boundary layer during each test day. The 
following subsections describe the equipment that was used 
for the experiments. 

2. sensors 

A complete set of meteorological equipment is used 
on the ship when it is outfitted for atmospheric research. 
The purpose is to obtain as complete a characterization of 
the atmospheric boundary layer (ABL) as is_ possible. 
Basically, there is need for a determination of the dynamic 
and mean parameters from the surface to the top of the ABL, 
often defined by a temperature inversion. For the purposes 
of the work described here, the dynamics of the layer are 
especially important because they are the driving mechanism 
behind dispersion from low altitude releases, and mean 
properties are needed in order to parameterize models in 
terms of readily measured quantities. 

A side view of the R/V Acania with the locations of 
the meteorological sensors is shown in Fig. 2-4. The ship 
has two masts, located on the bow, dedicated to the sensors. 
The foreward mast is on the tip of the bow and the sensors 
located there were at a height of 7 m above the mean water 
level; the second mast is 5 m behind the bow with sensors at 


20.5 m above the mean water level. This mast telescopes 
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Gown to a personnel platform so that the sensors can be made 
easily accessible. The platform also holds sensors that do 
not need to be elevated (aerosol counter and IR sensors). 

The ship is approximately 40 m long, 7 m wide, and 
only 7 m high (9 m at the ship's stack). The low profile 
and narrowness of the ship cause minimal disturbance to the 
air flow, making it ideal for overwater atmospheric 
research. The sensors on the high mast are well above any 
significant ship influence but there is some distortion of 
the flow at the elevation of the foreward mast. For this 
reason only data from the upper station are used in 
subsequent data analysis. Lower mast sensors are used as a 
backup in the event of an upper sensor failure. 

A summary of the monitoring equipment and associated 
meteorological parameters measured is given in Table 2-2. 
Details of the various pieces of equipment can be found ina 
previous report (Naval Postgraduate School, 1980). 

3. Data Acquisition Methods and Recording Procedures 

Four methods of data acquisition and recording were 
used: strip charts, analog tape recorders, computer 
controlled data acquisition and recording systems, and 
spectral analysis. 

Strip chart recording was used only for the acoustic 
sounder, relative wind direction and speed, anal the wind 
fluctuation signals from the hot films (TSI system). The 


internal strip chart was the only output available from the 
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TABLE 2-2 


METEOROLOGICAL MEASUREMENTS MADE ABOARD 
THE RV/ACANIA AND THE EQUIPMENT USED 


Measurement 
Relative Wind Speed 
Relative Wind Direction 


Air Temperature (T) 


Dew-Point Temperature (Tp) 


Sea-Surface Temperature (Tos) 


Wind-Speed Fluctuation (U') 


Three Axis Wind Velocity 
and Fluctuation (b) 


Ship Roll and Roll Rate 
Ship Pitch and Pitch Rate (b) 


Inversion Height (2};) 


Temperature Profile 


Ship Location 
Aerosol Content (a) 


Cloud Cover and 
Weather Conditions 
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Equipment 
MRI 1022 Wind System 


100 Ohm Rosemount. 
platinum resistor 
in a Gill aspirator 


General Eastern 1200 AP 
cooled mirror dew 
pointer, modified for 4 
wire resistance 


100 Ohm Rosemount 
platinum resistor and 
thermal ballast in 

a floating tube and 
Barnes PRT-5 infrared 
radiometer (a) 


TSI Constant Temperature 
Resistance Bridge and 60u 
platinum coated quartz 
resistance probes 


Kaijo-Denki Ultrasonic 
Anemometer 


Pendulums on the ship's 
pitch and roll axes 


Aerovironment 300 
acoustic sounder 


Radiosonde 


Loran C and Motorola 
Mini-Ranger III (b) 


Particle Measurement 
Systems Optical Counters 


Weather observations 


TABLE 2-2 (CONTINUED) 


Measurement Equipment 
Microturbulence Hot films at 60 feet 
Wave height Observation estimates 


a) Not used on BLM-4 
b) Used only on BLM-4 


acoustic sounder. The other strip chart data were seldom 
used for analysis. These recordings were made because they 
provide an immediate check on shipboard conditions. 

The analog tape recorder was essentially a back-up 
instrument. Every possible signal was recorded in this 
manner. If failure of the primary data acquisition 
equipment occurred, it was possible to retrieve data by 
using this recording. The temperatures were measured by 
resistors, which cannot be readily analog recorded. 

The central data acquisition components were the 
computer controlled data acquisition systems. Two were 
used: one dedicated to the ultrasonic anemometer and the 
ship motion sensors, the second devoted to obtaining 
meteorological data. A computer operated a scanner, 
voltmeter, and printer, and files data and calculated 
parameters on its internal cassette tape. 

The basic procedure for acquiring data for a given 
time period and averaging using a computer and scanner was 


straightforward and will not be described here. Only 
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average data and calculated parameters were stored to 
prevent using a large amount of computer memory and/or tape 
storage. All of the data and parameters were also printed 
at the end of an averaging period, providing a hard copy 
output and real time assessment of systems behavior. 

The actual averaging used was somewhat complex since 
both short term averages for turbulence parameters and long 
term averages for mean parameters needed to be obtained. 
Averaging periods used were 10 second and one-half hour 
intervals. A data acquisition cycle takes approximately 1 
second so that 10 readings were obtained for each short Ree 
average. All 10 second averages were held in computer 
memory until the end of the one-half hour period, when they 
and the mean data were averaged for the period. MThen both 
short and long-term averages were stored on tape and all 
long-term averages printed. 

True wind direction, corrected for ship's roll, and 
the true wind speed were obtained as short term averages 
from the meteorological data acquisition system. The 
ultrasonic anemometer outputs were processed to obtain short 
term averages of the three wind vectors, corrected for 
ship's pitch and roll. 

Spectral analysis had two functions: to determine 
the power spectral density of turbulence signals and to 


detect and identify system noise which would invalidate 
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results obtained by other acquisition methods. Normally, it 


was used on a regular basis only for the hot-film signals. 
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IiIt. DATA PROCESSING 


A. QUANTITIES OF INTEREST 

Data reduction is directed toward producing mean 
meteorological parameters and parameters which describe the 
turbulence in the marine layer. The mean quantities are 
easy to calculate from sensor response functions and the 
techniques need not be described here. Parameterizing the 
turbulence requires sophisticated techniques and is subject 
to error caused by measurement uncertainties and by 
misinterpretation of the measurements. In this section the 
several methods used to make these determinations are 
briefly described. The redundancy in methods provides cross 
checks on the results. 

Measurements made during the data collection expeditions 
of BLM 3 and 4 were made from a Single platform, so only the 
local turbulence was determined. Local turbulence is driven 
by a number of forces, which must be considered in 
developing any parameterization. Some obvious forces are: 
surface wind shear, convection, mesoscale activity, and 
swell and wind waves. It is important to establish which 
are the dominant factors to consider when developing the 


needed parameterization. 
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Specific quantities of interest can be categorized as 
either meteorological parameters or turbulence descriptors 
as listed below: 

Meteorological Parameters : 

True wind speed, U, 

Friction velocity, Ux, 

True wind direction, 6, 

Inversion height, 2j, 

Monin-Obukov Length, L = (T/kg) (Ux ?/H) , 

Surface Virtual Heat flux, H = PChUx Oy P 

Convective mixing velocity, Wx = [(g/T)H/z} 1/3 
Turbulence Descriptors : 

Standard deviation of longitudinal wind velocity, 9,,, 

Standard deviation of transverse wind velocity, 4,, 

Standard deviation of wind direction, Og 
In these equations, T = absolute temperature, g = accelera- 
tion due to gravity, k = von Karman's constant, 90 = air 
density, Cp = specific heat at constant pressure, Z; = the 
boundary layer depth, and 6yx*x = the virtual potential 
temperature scaling parameter. All of these quantities 
describe various aspects of conditions or turbulence in the 
atmosphere and are important for describing expected 
plume/puff properties. Turbulence generated by shear and 
buoyancy can be accounted for by the local meteorological 
parameters. Turbulence is also generated by mesoscale 


processes, which are difficult to account for with only 
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local parameters. The parameterization methodology is 
described in Chapter IV. 

The meteorological parameters were determined by local 
measurements of wind speed, air-sea temperature differences, 
Gew-point temperature, and inversion height. The 
Monin-Obukhov length characterized the state of the surface 
layer based on the hydrostatic stability. The surface layer 
momentum, heat, and water vapor fluxes are needed as a 
measure of forcing in the boundary layer. The scaling 
parameters can be used to calculate the surface layer fluxes 
of momentum, heat and water vapor. In a well mixed boundary 
layer, Wx is the appropriate velocity to use to determine 
the rate of mixing throughout the convective boundary layer. 
The stability, friction velocity, dissipation rate, and 
Mixing velocity were all determined using the 
bulk-aerodynamic method. The method is valuable because it 
makes use of readily measured quantities. Specifically, it 
is based on differences in mean temperature, wind, and water 
vapor at the surface and at some reference height. For a 
more complete description of this method, see Schacher, et 
al. (1986). 

The next section describes the quality control and error 
checking procedures performed prior to further data 
manipulation. This is followed by a description of the 


process used to derive standard deviations of horizontal 
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wind velocity components in the downwind (u) and crosswind 


(v) directions. 


B. WIND DATA CALCULATIONS 

The objective here was to process the BLM data base in 
order to determine standard deviations of horizontal wind 
velocity components in both the longitudinal and transverse 
relative wind directions (94 y). This information was then 
plotted and used as a data base to which a theoretical 
parameterization was empirically fit. In Schacher et al. 
(1986) a parameterization. for overwater wind direction 
variability was developed by using the standard deviations 
of wind direction (dg). The vector method used here to 
derive standard deviations differs from the method used by 
Schacher et al. (1986) since it conceptually treats 
contaminant motion differently. The two methods are 
compared later in this report. 

During each step of the data processing phase, 
provisions were made to double check the resulting values 
Since significant errors are easily made processing 
turbulence data. Only the BLM 3 and BLM 4 See collected 
0, data, therefore, only the data from those two cruises 
were used here. The deletion of BLM 1 and BLM 2 cruise data 
reduced the number of data sets from 859 to 705. The 
remaining data sets were quality checked for validity. In 
each data set there were between 39 and 150 wind component 


data elements. There was a possibility that some data were 
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lost when transferred to the currently used format. This 
was confirmed by printing out the number of data elements in 
each set and comparing them to the archived records of what 
each data set should contain. Data sets with too few wind 
components were not used. In the final quality control of 
the BLM 3 data set, 30 out of 401 sets were discarded. This 
accounted for 7.5% of the sets. In the BLM 4 data, 16 out 
of 304 sets, 5.3%, were discarded. Once a "clean" data base 
of 659 data sets was established, test processing of one 
data set was done. This same test data set was processed by 
hand to verify the methodology. After this successfully 
verified the technique, all data sets were processed. 

Puff motion may be described by transport and diffusion 
characteristics. Transport can refer to either the mean 
trajectory of the puff over a time period, or the meander of 
the puff center of mass about the mean trajectory within the 
wind field. Relative diffusion refers to the puff growth 
about the puff center of mass. These puff mechanisms are 
influenced by differing scales of turbulence, where 
turbulence refers to any variation of wind flow from the 
mean. Large turbulence scales, much larger than the puff 
size, move the whole puff as an intact entity in the wind 
field, resulting in meander. Small turbulence scales cause 
relative motions of parcels of air that are internal to the 
puff. This relative diffusion causes the puff to grow about 


its center of mass. 
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Most data sets contained enough data elements’ to 
determine an average particle trajectory over a 30-minute 
period. Of course, using only the 30-minute average does 
not fully describe what would actually happen to a small 
puff release. As described above, meander and diffusion 
affect a puff as it flows downwind. When a puff is initially 
released, the size of the puff is small. As the puff moves 
downstream within the wind field, the puff changes due to 
forces internal to the puff, or due to external wind shifts. 
If a puff of length '‘l' moves distance 'd' downwind in time 
't', then puff growth during small length scales take place, 
where 1 = qd. Large length scales refer to puff growth that 
occurs where l << d. These changes can occur within 
different time scales. Small internal changes may occur at 
short time intervals within an area relative to only the 
puff's instantaneous location in space. Large external 
influences may occur over longer time intervals, such as a 
Wind shift, which causes the whole puff to change its actual 
location relative to the ground. 

In order to study puff growth due to wind variability, 
the BLM data was processed into specific time averages, then 
analyzed. Specifically, the wind speed and direction infor- 
mation were vectorized and processed to provide data element 
sets containing U, 6, g, Sy for 1-, 3-, 10- and 30- minute 
time averaged intervals over each 30-minute period 


covered by a data set. The vectors were computed using 
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Yi = Uy * cos(9;), 
Xi; = Uj * sin(%;), 


and then summed into time bins. (Recall that the original 
data contained wind speed and direction.) The wind vector 
components were yj; for the north, and xj; for the east 
components, respectively. Note that the angular origin is 
along the cartesian +Y axis to correspond with true north 
(000°). 

The shipboard wind sensors computed time averaged data 
elements at about 14-second intervals. For the purpose of 
this research, each one-minute time interval contained four 
data elements yielding 56 -second "minutes." Each data 
element was summed into the appropriate time averaged bin 
such that there were 4 elements/1 minute; 12 elements/3 
Minutes; 40 elements/10 minutes; 120 elements/30 minutes. 
For each time bin which contained the appropriate number of 
data elements, the average x, y vectors were computed. The 
average X, Yy vectors were converted back into an average U, 
6 for the time bin. Standard ey eiions from each bin were 


calculated using the formulas: 


ieee acim) ian. 
E N-1 
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and 


oe [Fepe | 
V N-L P 
where 
uj; = Uy * cos(% - Genin, 
and 
Vi = Ui * sin(6; - 6(t)) , 


where a(t) was the average wind direction for each of the 
four appropriate time bins. 

Counters were used to subtotal the number of times each 
bin was accessed so that after the a were processed, the 
mean wind speed (U), direction (co), and standard deviations 
(Su,v(t)) were computed for each of the time periods. The 
values for the four time averages were then taken and made 
into grand averages of 30 minutes. For example, 30 one- 
minute wind speed average values were determined from one 30- 
minute data set, and all 30 were averaged to get the overall 
11-minute average for the entire 30-minute data period. Due 
to arithmetic properties (e.g., commutative, associative, 
additive), the 30-minute grand average wind speeds and 


directions for all four time averages were very similar. 


3:2 


The oy,y(t), U,@ values were read into an output file for 
comparison with the U, 8, dy, and standard deviation of the 
wind direction (op) values generated by Schacher et al. 
(1986) and for further processing. 

A major strength of the vector method was that it 
provided a better physical description of wind flow in the 
boundary layer. This method described an average particle 
trajectory path based on averaging a series of vectors. The 
vectors accounted for both direction and magnitude of a 
particle's velocity, and therefore more accurately described 
the physical processes occurring as the particle moved. 
There was a high confidence level in the accuracy of 
information generated by the vectorization method. The 
scalar method used previously was’ strictly a scalar 
averaging of U,@ without regard to the occasional large 
variations in the vector magnitude and direction within a 
data set. By only taking either angular averages, or only 
speed averages, the scalar method ignores the effects of 
interaction between sequential vectors having large 
differences in either/both direction or magnitude. The 
result was that an average particle trajectory path was less 
accurately described than by the vector method, making the 
previous results less usable for our purposes. Fig. 3-1 
illustrates an example of each method's averaging of the 
same two vectors. As can be seen, the vector method 


produces a more accurate average than the scalar method. 
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a= 76 
p =1 28° 
o= 772° 





Vector method utilizes vector addition, representing 
true transport. 


Scalar method utilizes straight, unweighted averages 


c= 6 
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Y= 95° 
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Fig. 3-1 Comparison of Vector versus Scalar 
Method Used to Average BLM Data. 
Part a Shows the Results of Vector 
Averaging, Part b Compares the Two 
Methods. 
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The scalar method will produce the correct U if the angles 
are the same and the correct\theta if the speeds are the 
same, but makes incorrect assumptions in the relationship 
between the two vectors when this is not the case. Example 
3-1 lists U,® for a data set showing the differences between 


scalar & vector averaging for a 30 min time average. 


Example 3-1 


Nu= 121 Total= 121 


a 2 

NPUT I J B(I,J) 

tas 1827 2m 2032 3 4 Bea 4 4 oe 
2 202 22 384 a 714 4 2 1069 
; 8 435 2 3 438 3 3 439 4 3 435 
1 @2 AS 2 4 28 3 4 29 4 4 75 
UTPUT IJ B(I,J) 

; a 20 a 20 ot 20 4 | 20 
ft 2 16 2 2 23 3 2 S52 4 2 79 
Ps 427 2 3 427 3 3 424 4 3 420 
| 17 ee 24 3 4 4 4 73 
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Once processed, all clean data sets (659) were compared 
to the corresponding Schacher (1986) data sets for 
comparison of each's results. In virtually every case, the 
average wind speeds were within +0.2 m/sec of the earlier 


values. For any given data set processed by the vector 


“method, the average wind direction value for each of the 


four time averages was within several degrees of each other. 
This was not the case for data sets processed with the 
scalar method. Differences up to 180° occurred between the 


four time averaged values in some data sets. In comparing 
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the results from both methods, the wind directions were 
rarely identical, but differences between the 30-minute 
averages were <+20° in 98.5% of the cases. Differences 
between the il-minute time averages were <+20° in 90.9% of 
the cases. Example 3-1 shows a typical case where there was 
a large disparity between the two methods in the 1-minute 
averages, but showed near agreement for the 30-minute 
average. 

A possible reason for the large fluctuations in average 
wind direction by the scalar method was not making proper 
corrections during computations where the wind directions in 
a data set were scattered around north (000°). An overly 
simple example is 6 = Hesse + 005°)/2 = 180°, which is 180° 
away from the correct average of 000°. In 76.6% of the data 
sets where the scalar method average wind direction value 
— greater than 20° different from the vector method's 
average wind direction value, the wind data were in the 
quadrant centered around north (000%). This provides a 
strong indication that incorrect adjustments were made for 
wind direction values centered around north. Example 3-2 
lists the wind directions in a sample data set and provides 
the scalar wind direction averages for each of the four time 
averages used. As can be seen in the example, the 1- and 3- 

minute time averaged values are very close to those derived 


by Schacher et al. (1986). The vector method 
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a 


automatically took this possibility into account to ensure 


these errors did not occur. 


Example 3-2 
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Other data comparisons included converting Oy based on 


the vector method into g, values, using some of the many 


8 
data sets which contained U, 4, oa, values that were in close 
agreement from both methods. This was done using oy = 
arctan(U * Tp )s where o 6 is in radians. Numerous random 
values of oy, were converted into og and compared to the 
Schacher et. al. (1986) values. In all cases tested, the 
two oO, values were in close enough agreement to indicate 
that the two methodologies produce essentially the same 
results. 


Lastly, a comparison of o,, values was made. In all 


cases examined, the scalar method's 0, values and vectorial 
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method's oy,y values were consistently a factor of two 
times within each other. As expected, the O's increased 
as the time average increased (1.e.: 94 < 0% < Og < 630) - 
The rate of increase for the o's was about the same by 
both methods. Based on hand calculations, and comparisons 
to the Schacher et al. (1986) data base, it is concluded 
that the vectorized data are valid and useful. The next 
chapter applies turbulence theory to the development of a 


horizontal wind variability parameterization. 
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IV. TURBULENCE THEORY 


A majority of the theory used here to develop a 
parameterization is taken from Schacher et al. (1986) and 
modified to meet the specific needs of this work. The 
following treatment follows from the results developed by 
Hojstrup (1982), where expressions were developed for the 
turbulence intensity as functions of uy and wy. Shear and 
buoyancy production were considered as the two sources of 
turbulence and were expressed as functions of the convective 
mixing velocity (wx) and friction velocity (ux). Their 
contributions to the velocity spectra were modeled and the 
spectra integrated over the appropriate frequency range to 
determine the velocity variances. The following is a brief 
description of how the Hojstrup treatment is used to develop 
algorithms for parameterizing overwater horizontal wind 
variability from our data. 

Scaling parameters were developed as a preliminary step 
in developing the equations. For some quantity xX, the 
relation between the value at height z (Xz) and surface 


values of X, and its scaling parameters is 


Xe - Xe = (Xx/Oyk) (ln (2/Zoy -Vy(L)]. (5) 


ao 


where the von Karman's constant, k, is 0.35. Rather than 
use the more conventional k = 0.4, k = 0.35 is used to be 
consistent with the calculations of Ux and wy from the data. 
Zo is the roughness length and y the stability correction 
function to the logarithmic profile, which is normally 
written as a function of the Monin-Obukov length (L), as 
indicated. The turbulence diffusivity ratio is a. 

The horizontal, transverse velocity spectrum has two 


components 
£S(f) = A(B,)we? + B(B5)Ux2, (6) 


where f is the frequency, S(f) the spectral intensity, and 


A(Bz) and B(B 9) are functions of the variables 
By = £(2;/U), Pipe tay ee, (7) 


where Z; is the height of the lowest inversion, and Z is the 
measurement height. The functions, A and B, were determined 
by Hojstrup by matching model results to the Kansas (Kaimal 
1978) and Minnesota (Kaimal et al. 1976) data, which gives 
spectral intensity for each component A, and B. A complete 
description in the methodology is contained in Hojstrup 
(1981). 

Upon integration of the spectral intensity, the velocity 


variances were found by Hojstrup to be : 
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Q 
Nw 


ye = 0.6k2/3 wa? + 4.8(1-2/2;) 2U42/(1415Z/24) 72/3 , = (8) 


0.7k2/3wy2 + 2.7(1-Z/Z4) 2U42/(1+2.82/Z;) 72/3 . (8a) 


Q 
< 
N 
| 


The first terms are normally written as 
0.6(-Z4/L) 2/3U,y? , (9) 
and 
0.7(-Z;/L)2/3u,2 . (9a) 


These equations are applicable only during unstable 
conditions, since that is when buoyancy driven convection 
occurs and wy is defined. During stable conditions, ws = 0, 
and the first terms in Eqs. (8) and (8a) are zero. 

It is apparent that the first term in Equations 8 and 8a 
contains Z;, through wy, because the strength of convective 
mixing depends on the mixing depth. The second term 
contains Z;, through the ratio 2/Z;, which is present to 
correct for the decrease in surface shear produced 
turbulence with height above the surface. The correction 
factor is 1 at the surface, since the numerator in the 


second term of Eqs. (8) and (8a) approaches zero there. 
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Since measurements are all within the surface layer, this 
factor should not be appropriate, therefore, the height 
correction in the second term of Eqs. (8) and (8a) are 
dropped. However, for a very low inversion it may play a 
role in indicating suppression of turbulence at the 
measurement height (20 m). Also, it is recognized that the 
overland spectra from which Hojstrup's model were obtained 
IMay not adequately represent the overwater case dealt with 
here. Thus, the numerical factors are incorporated into the 
production mechanism parameters C, and Cy, and the 
correction factor in the second term is set equal to 1.0. 


The resulting equations are 


Q 
li 


or 


Oy vy SiG >247 E)-? eee (iy 
where C, and C, are constants to be determined from the 
overwater data. Note that the correction factor in the 
second term was suppressed; if necessary it could be added 
later. For simplicity C,' = Cyk2/3 will be used. 

The assumption is made that turbulence behaves in 
Similar patterns over water and over land. Similar 


production mechanisms are expected to be valid for both 
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regimes, although the magnitudes may be different. bt 
spectral intensities were plotted for both regimes they 
would be expected to have similar shapes. Both regimes need 
similar algorithms to describe the geophysical processes 
contributing to turbulence and as such, both depend on 
scaling parameters. Hojstrup (1982) has already shown this 
for the overland regime. It is a goal of this work to 
extend the methodology of the previous work to the overwater 
regime. 

It may be convenient to write the scaling velocity as a 


function of the wind speed and drag coefficient, 
U. = cl/2y (13) 


Following Eq. (5), the drag coefficient can be written as 


cl/2 = cy[1 - Cyt 2 ¥/k]~7 } (14) 
with the neutral stability drag coefficient 

Gy / ¢ ="k/1n (2/25): (15) 
For unstable conditions (Businger, 1973) 


v= 2 In{(1+x)/2] + 1In({(1+x2)/2] - 2 tan71x + pi/2 , (16) 
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with 
x = (1 = 152/L)+/4 
and for stable conditions 
W = -4.7(Z/L) . (16a) 


At this point, the parameterization reduces to using 
determined values of the stability (Z/L), inversion height 
(Z;), wind speed (U), and roughness length (Z,) or neutral 
drag coefficient (Cy) and finding the values of C, and C,,. 

The roughness length, or the drag coefficient, over 
water depends on the wind speed through wind-wave 
interaction. The Garratt (1977) formulation for the neutral 


stability drag coefficient, 
103Ccy = 0.75 + 0.067U , (17) 
was applied to the buoyancy production component of the 
regression equation, so that the buoyancy production is 
represented as. 
Us? = (7.5 x 1074 + 6.7 x 107°U)F(L)UZ , (18) 


with 
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F(L) = [1 - Cyl/2 0 /kj}72 . (19) 


This is not a closed form solution since L depends on U. 
All coefficients are determined by an iterative calculation 
when the data are obtained. 

Rather than do another iterative calculation for the 
results presented here, and in order to clarify the results, 
F(L) was calculated using the old value of L and Equation 
HO. Note that yw for stable conditions was given by -4.7 
Z/L. At wind speeds of 10 m/sec or higher, Z/L is near zero 
Since the air-sea temperature difference never can be large 
enough for buoyancy to overcome shear production at such 
high wind speeds. 


Eqs. (11) and (12) can now be rewritten as 
Sy = Cy Wa + CyP(L) (7.5 x 10-4 + 6.7 x 10°7°U)U2, (20) 
where 
Ga = Fee 90. 497°C, 
and, 
eee 2a 0) -/ eer (Ly (75 x 10=- + 6-7 x 107°U)Ue. 


(21) 
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The first form was used for wy and U parameterizations and 
the second for stability parameterization and for the 
effects due to inversion height. Recall that the inversion 
height correction to C, (see Eq. (8)) may be needed for low 
inversions. Another form of Eq. (21) is useful for 
examining the dependence on stability. It was found by 


writing the term in square brackets as 
Cyl (24/2) (-Z/L) 273 + c, . (21a) 


In Eqs. (21) and (21a) both L and U are left as parameters. 
Stability did not appear as a sufficient variable. 

It is reemphasized that none of the above addressed 
sources of turbulence other than shear and buoyancy 
product jon. The mesoscale process could not be 
parameterized solely with C,, and C,, so another term needed 
to be added to Eq. (20). Mesoscale production is a function 
of differential heating between the land and sea which 
contributes to low wind speed turbulence. The contribution 
due to thermal convection decreases with increasing wind 


speed such that 
Su,v (mesoscale) = yun : (22) 
If N = 2, Cys can be interpreted as the square of an 
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Mixing velocity, ws, and the friction velocity U,. 
Understanding this new velocity is not easy, and here only 
the new constants were evaluated with no explanation of 
their meaning. The final algorithm used for the 


parameterization became 


Gy vy? = 0.497Cywe? + CyF(L) (7.5 x 1074 + 6.7 x 1075u)U 


+ Cyo/UN , (23) 


where the three fitting parameters were functions of the 
averaging time used to evaluate the variances. The first 
term represented contributions due to shear, the second term 
accounts for buoyancy, and the third for mesoscale 
production mechanisms. Stability was taken into account 
theough F(L). Based on Schacher et al. (1986), it was 
determined that oversea conditions did not support stability 
conditions far from neutral for high winds, and that the 


stability correction was not needed, therefore 
Fiije= @h0 . 
The parameterization described in Chapter VI was 


performed using this algorithm to obtain an empirical fit of 


the horizontal wind velocity component variances to the data. 
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V. DATA PROCESSING RESULTS 


A. GENERAL 

Parameterization of the wind variability required 
segmenting the data into wind regimes. The situations were 
different when the wind was dominated by synoptically driven 
northwest flow than when it was dominated by the sea-breeze 
cycle. As pointed out in Chapter III, buoyancy, shear and 
mesoscale production mechanisms are different for various 
meteorological conditions, which affects the degree to which 
wind variability exists under a given set of environmental 
conditions. 

In Schacher et al. (1986) multiple data subsets were 
generated for well known environmental conditions such that 
a given forcing was dominant and a single parameter could be 
identified as the dominant factor in parameterizing the wind 
variability. One subset was the stationary wind regime 
where the one-hour averaged wind direction remained constant 
within about 40 degrees and the wind speed was fairly 
constant, for which mesoscale activity was minimal. In some 
cases clean data subsets were not available. When more than 
one parameter was the dominant factor, self-consistency 
among the dependences was examined to determine dependences 
on the individual parameters. The analyses focused on 


dependences on the wind speed, convective mixing velocity, 
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and stability, with inversion height and wind direction as 
important auxiliary parameters. The dominant parameters 
considered for parameterizations were established after 
examining the following variables: 

averaging time, 

onshore/offshore flow, 

stability, 

Wx, 


Zi- 


B. SCATTERPLOTS 

The initial set of plots were scatterplots for Ou,v/U 
versus wind speed (U). The reason for dividing the standard 
deviations by wind speed was to non-dimensionalize the 
values and to compare with later theory. The scatterplots 
were produced for both BLM-3 and BLIM-4 data, for each of the 
four time averaged periods of 1, 3, 10 and 30 minutes. A 
full set of scatterplots is contained in Appendix A. The 
large variability that was normally present in fluctuation 
data is evident. The standard deviation values increase as 
the averaging time increases. The 1- and 3-minute time 
averages have a oy,y/U value less than 0.3, increasing to 
0.6 for the 10-minute average, and to 1.2 for the 30-minute 
average. This is a four fold increase in the amount of wind 
speed variability just from using longer averaging times for 
the same data. As the time averaging period increased, more 


time elapsed within which to collect wind direction data. 
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This allowed dramatic shifts in wind direction to occur 
within a data set, which in turn, resulted in larger 
deviations from the mean wind. In a time period of 30 
minutes, it was possible for the wind direction to shift up 
to 180° for some conditions. Over a short time period of 1 
minute the amount of wind shift was much less, hence, the 
smaller standard deviations. 

For all four time averages the oO, variability is less 
than the corresponding 09, variability. Even though there 
is considerable scatter, trends are evident in the plots to 
demonstrate that there is consistency in the data collected 
@uring both cruises. In what follows, the term 
"variability" will be used to mean the wind speed standard 
deviation in either the longitudinal (downwind) or 
transverse (crosswind) directions. It was difficult to 
obtain quantitative information from scatter plots. They are 
mainly useful for verifying consistency and trends within 
the data. In the next section plots of the averaged data 


are discussed. 


C. AVERAGED DATA PLOTS 

The data used in the scatterplots were averaged over 
discrete wind speed ranges of 1 m/sec increments from 0 to 8 
m/sec; 2 m/sec increments from 8 to 12 m/sec; and a 3 m/sec 
increment from 12 to 15 m/sec. The averaged data plots for 
unrestricted conditions are shown in Appendix B. These 


figures also contain a print out of the results averaged 
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over all data, giving the value of the variability for the 
center of that range; the number of points found in that 
range; the mean standard deviation of the horizontal wind 
speed over the averaging period indicated (sigma); and the 
standard deviations of the data about the mean sigmas. Ina 
following section an algorithm is developed which is 
empirically fitted to these averaged plots for unrestricted 


conditions. 


D. PRODUCTION MECHANISMS 

In the interest of minimizing redundancy of information 
between here and Schacher et al. (1986), only a_e short 
overview is presented of the earlier conclusions regarding 
dominant meteorological parameters used _ to develop 
parameterizations for wind direction standard deviations. 
Those findings were used here as a stepping off point, 
forming the basis for the approach used to develop 
parameterizations of horizontal wind speed variability. 

Criteria for stationary conditions were: the one-hour 
average wind direction remained within a 40 degree sector 
and the wind speed was fairly constant. Such conditions 
only occurred when the wind was fairly strong and from a 
westerly to northwesterly direction. Mesoscale forcing was 
largely absent during stationary conditions since the strong 
steady wind flow negated any mesoscale effects due to 
thermal convection. This factor allowed for the buoyancy 


contribution to be determined. 


oe 


Stability was a function of the air-sea temperature 
difference. Due to the small temperature differentials, 
only a small range of stabilities was encountered over 
water. Since buoyancy production occurs only when there is 
upward heat flux, it was absent during stable conditions. 

The shear production contribution was isolated during 
stationary, stable conditions and was found to be dominant 
at wind speeds above 7 m/sec. Below that speed, buoyancy 
production was apparent. 

Based on Schacher et al. (1986) no apparent difference 
in the turbulence was seen for onshore and offshore 
conditions. Convection normally played a minor role in 
turbulence production aS waS evidenced by the turbulence 
being slightly greater during unstable conditions. This 
weak dependence on stability led to utilizing stationary 
conditions to examine turbulence dependence on wy. Large 
values of variability at low wind speeds were associated 
with non-stationary conditions during transition periods in 
the land-sea-breeze cycle. The forcing for the cycle was 
differential heating between the land and the sea; a 
mesoscale process. Low wind speed turbulence produced other 
than by buoyancy are referred to as "mesoscale" production. 

The BLM data were separately processed and averaged for 
stationary and for stable conditions so that the effects of 
the various production mechanisms on wind variability could 


be individually isolated. Appendix C contains averaged plots 
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of oy y/U versus U for stable conditions for all four time 
averaged periods. Stability calculations included the shear 
and buoyancy production mechanisms, which were a function of 
mixing depth (Zj;). Since buoyancy production was minimal 
during stable conditions, the plots in Appendix C for stable 
conditions give indications of shear production 
eontribution. It can be seen from these plots that the 
variability decreases as wind speed increases in all cases, 
as it does during unrestricted conditions. Also, the values 
of Oy,v/U increase from 0.15 to 0.75 for U < 1.0 m/s, as 
the averaging time increases from 1 to 30 minutes; a 
five-fold increase. Appendix D contains averaged plots of 

%,v/U versus U for stationary conditions for all four time 
averaged periods. As previously mentioned, this situation 
allows for an estimation of buoyancy production since shear 
and mesoscale production are minimal. In all eight cases, 
the plots show that the contribution is relatively constant 
with respect to wind speed, and that the Sy,v/U values are 
approximately the same. 

Schacher et al. (1986) determined that although 
stability was a good parameter for shear and buoyancy 
production, it did not account for mesoscale production. 
Based on the previous work, the optimum conditions for 
parameterizing mesoscale production was stable conditions, 


at low wind speeds (< 7 m/sec). 


23 


Chapter VI discusses how the theoretical algorithm is 
used to determine regression variables and generate fitting 
curves for the averaged data plots. Conclusions are 
provided with implications for future applications of the 


parameterization results. 
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VI. PARAMETERIZATION RESULTS AND CONCLUSIONS 


A. PARAMETERIZATION RESULTS 

To parameterize horizontal wind velocity variability, 
wind data were decomposed into longitudinal (downwind) and 
transverse (crosswind) velocity components over different 
time averages. Standard deviations of horizontal velocity 
components were then computed and plotted. Horizontal wind 


variances were then parameterized using the regression 


equation 


Cro] OMOTC wa + C,(7.5 x to“ 


+ 6.7 x 107°U)U2 + c,./UN , (23) 


where the regression variables representing shear (C,), 
buoyancy (C,,), and mesoscale (C,.,) production are different 
for the u, and v_ directions. The three regression 
parameters are functions of the averaging time used to 
evaluate the variance. 

Presentations of the data in previous sections of this 
thesis and by Schacher et al. (1986) allowed some 
conclusions to be drawn about the mechanisms driving 
horizontal wind variability and the methodologies needed to 
determine the correct parameterization. Three production 


mechanisms were identified from the data: namely buoyancy, 
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shear and a larger scale forcing associated with mesoscale 
processes. The previous work showed that it was not 
possible to unambiguously separate the mechanisms. Certain 
associations between these production mechanisms and 
associated conditions were used to separate the processes 
for purposes of developing the parameterization. 

Shear production dominated at high wind speeds (> 7 
m/sec) and was insignificant at low wind speeds. The shear 
production contribution was determined from stationary and 
stable thermal conditions at high wind speeds. Mesoscale 
forcing was absent during stationary conditions and buoyancy 
production was absent during stable conditions so that shear 
production could be isolated and individually parameterized 
using the second term of Eq. (23). Based on the previous 
work, it was determined that oversea conditions were near 
neutral for high winds. Thus, it was assumed that the 
stability correction was not needed for the values of Cy, 
Cy, Cmsg to correctly fit the data. The result was to use 
F(L) = 1.0 in Eq. (23). 

Buoyancy production was apparent at low wind speeds (2g 
m/sec) and decreased inversely proportional to wind speed 
(1/U). Good results were expected when using stationary 
conditions at low wind speeds to determine buoyancy 
contributions. This was due to minimal mesoscale production 
for stationary and minimal shear production at low wind 


speeds. Figs. 6-1a,b show plots of o,/U vs. U for the 
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1l-minute average during unrestricted and stationary 
conditions, respectively. When comparing the two figures, 
it is seen that by isolating stationary conditions, buoyancy 
effects at low wind speeds and shear effects at high wind 
speeds may be isolated, and thus quantified. 

Mesoscale production was deduced by subtracting the 
shear production contribution from non-stationary, stable 
conditions. The mesoscale processes dominated at low wind 
speeds (<5 m/s) under these conditions, since bouyancy and 
shear production were minimal. 

The parameterization for horizontal wind variability was 
developed by fitting Eq. (23) to the data, then adjusting 
Cy, Cyr, Cmg and N to find the best fit. This was done for 
all four averaging times of 1-, 3-, 10-, and 30-minutes for 
both 9,, and q,. 

Appendix E contains the empirical fits of Eq. (23) to 
the o, y/U vs. U data for stable, Stationary and unrestrict- 
ed conditions. Sorting on wy ranges of 0-0.2, 0.2-0.4 and > 
0.4 m/sec was done for the latter two cases and the values 
of the fitting parameters used are shown on the graphs. 
Figs. 6-2a,b, and c provide examples of wy sorting for 9%,,/U 
vs. U during stable, stationary and unrestricted condition, 
respectively. The ranges of 0.0-0.2, 0.2-0.4 and >0.4 m/sec 
are represented as 0, 2 and 4, respectively. These 
particular conditions were chosen because a) ws = O for 


stable conditions and Cns can be determined, b) 
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mesoscale production was low for stationary conditions and 
Cy can be determined, and c) the whole procedure can be 
checked for unrestricted conditions. For each set of curves 
the parameters were evaluated by the following procedures. 
Figs. 6-2a-d contain plots of o,/U vs. U for the l1-minute 
time average, demonstrating the sequence of these 
procedures. 
1. Choose an initial value of C,, based on the difference 
between the values for stationary, and _ stable 


conditions at high wind speeds. 


2. Adjust Cypsy and N to fit the low wind speed portion of 
the stable data, letting C, = 0. 


3. Adjust C,, to correctly fit the high wind speed regime 
during stable conditions, letting Cy = 0. Fig. 6-2a 
contains an example of the adjustments performed in 
steps 2 and 3. 

4. Choose C, and adjust to fit the high wind speed regime 
during stationary conditions, letting Cyo = 0. Fig. 
6-2b contains an example of this step. 

5. Check the results by checking the fit that all 
parameters combined give to the unrestricted data. 
Fig. 6-2d contains an example of the end results for 
this particular case. 

There are two fitting curves for those cases where C, <> 
O. The lower curve is for wy = 0.1 and the upper curve is 
for Wx = 0.5 m/sec. The fitting curves turn upward below 2 
m/s as the mesoscale contribution (Cyo/UN) dominates over 
the shear and buoyancy contributions. Fig. 6-3 shows the 
0.,/U vs. U plot for the 30-minute time average (the values 
representing o,,/u for the 30-minute time average will 


hereafter be referred to as 0439). Using the values of the 


regression variables presented in Appendix E, o 307 may be 
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calculated using Equation 23, assuming wy = 0.1 m/sec, such 


that 


Suso’ = 08497(0-3)0.1)-¢ + (0.01) (7.5xH0me 


+ 6.7x107°U)U2 + 0.28/U2 . (23a) 


Applying Equation 23a during various wind speeds, the values 
of the three subcomponents representing shear, buoyancy and 
mesoscale production mechanism contributions are calculated 
and displayed in Table 6-1 along with the o 439% and oy3o 


values. 


TABLE 6-1 


COMPARISON OF PRODUCTION MECHANISM CONTRIBUTIONS 
BASED ON EQUATION 23 WITH RESULTANT Oo, AT 
VARIOUS WIND SPEEDS 


U (m/sec) Shear Buoyancy Mesoscale 51307 130 
0.1 IR a SE 7.5E-8 28.0 28.0015 S23 
On Re: LE =6 Ieee 1.1215 1.06 
5.0 1h 2.7E~4 1.3E-2 1.477E-2 0.114 

11.0 ao Lae 


Teo o—s 2s J30e35 Sy 1S 7. SEs 


Table 6-1 shows that as expected, at wind speeds below 5 
m/sec, mesoscale production dominates. 

The values of the empirically derived regression 
variables are given in Table 6-2. For the purpose of 
comparison, the values determined by Schacher et al. (1986) 


Lor oe are included. The shear and buoyancy values in 


se 
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Table 6-2 are about the same size for both %, and 4. This 
result is expected based on previous discussions by Panofsky 
and Dutton (1984) and Hojstrup (1982). No direct comparison 
should be made between the magnitudes of the regression 
variables since the actual significance of this work lies in 
analyzing the values of On v2. The values for N are similar 
for 9°, but smaller for °%,. 

It was found that the best value of N for the mesoscale 
term was dependent upon which wind component was examined. 
In the transverse direction (v) N = 2 for all four time 
averages. This was expected since Schacher et al. (1986) 
determined N = 2 for 9, and Sy, is directly related to %. 
Fig. 6-4a shows the o,/U vs. U plot for the 30 minute time 
average. Fig. 6-4b shows the corresponding SC, vs. U plot 
for the 30-minute average from Schacher et al. (1986). 
Table 6-3 shows sample values picked off of the graphs. The 
"0, conversion" column contains the equivalent values of 0, 
converted too, using the formula 


6 


{tan[(9 ,/U)*Uj}/U = O, (radians) ; 


and compared to the corresponding Op values from Schacher et 
al. (1986). As seen in Table 6-3, the converted oO, values 
are very close to those produced by the previous work, 


indicating consistency between the two methodologies. 
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Fig. 6-4a 06,/V Plot for 30 Minute Time Average, 
Unrestricted Conditions, Using BLM 
3 and 4 Data 
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Oo, vs. U Plot for 30 Minute Time Average, 
Unrestricted Conditions, Using BLM 1, 2, 
3 and 4 Data (from Schacher et al., 1986) 
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TABLE 6-3 


COMPARISON OF og VALUES CONVERTED FROM 
oy AND FROM SCHACHER ET AL. (1986) 


U (m/sec) G/U On (Schacher) Og (conversion) 
Ih 0.34 Z1.3° 2s 

ce 5 0.09 Eee Sao 
11.0 0.055 23 Be 3.6° 


In the longitudinal direction (u), many values of N were 
tried, with N = 1.33 yielding the best fit for the 1, 3, and 
10 minute time averages. Figures 6-5a,b compare the 
fitting curves based on N = 2.0 and N = 1.33, respectively. 
It is seen that the more gradual slope of N = 1.33 is 
required to properly fit the low wind speed data. Oddly 
though, N = 2 provided the best fit for the 30-minute time 
average. Since previous works (Schacher et al., 1986; 
Hojstrup, 1984; Panofsky and Dutton, 1984) utilize 30-minute 


time averages, it is necessary to emphasize findings based 


on the 30-minute averaging which shows N = 2.0 for both the 
u and v directions. The implication is that mesoscale 
effects are about the same in both directions. The 


different values for each fitting parameter suggest that the 
production mechanisms affect puff growth differently in the 


downwind than in the crosswind directions. 


69 


Sigma-U 7 U 


Cu = .@2 1 min Ave 
Cu = .00? | Total 
Cms = .@25 | 

* MESO SLOPE = 2 


Sigma-U 7 U 





Cw = .@618 l min Ave 
Cu = .4G?. Total 
Cms = .811 

MESO SLOPE = 1.33 





Q | 5 12 iS 
/ WIND SPEED (m/sec) 


Fig. 6-5 Comparison of Empirical Curves Using 
Mesoslopes of 2.0 Versus 1.33 
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B. RATIOS 


Panofsky and Dutton (1984) proposed the _ simple 


expression 

QA * Ux , (24a) 
and 

OF = B* Us , (24b) 
for purely shear driven, mechanical turbulence. A and B 


were presumably constants for a given surface roughness, 
assuming convection was unimportant. He also assumed that 

Oy was slightly larger than On. Based on various 
experiments over flat terrain the mean values for the 
constants A and B, a 30 minute time averages were 


A= 2.39 + 


JL 
© 
Oo 
UW 


and 
be-weI2 - 0.05 , such that, 
hy Bemeda25 + 0.05... 
These values agreed well with corresponding wind tunnel 
measurements. The assumption here is that turbulence is 
shear driven and is quantified by the oO/uy, ratios, 
therefore, only uy, is a function of the generation 
mechanism. Under this assumption, A and B_ should be 
universal constants based on the integral of the velocity 


spectra over the time scale of interest. 
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Based on observations at various locations, it was seen 
that as the surface roughness increased, the values of A and 
B increased. Conversely, over open ocean having a small 
surface roughness, it would be expected that A and B 
decrease, perhaps being several times smaller than the above 
values. Table 6-4 shows the values of A and B based upon 


the BLM data. 


TABLE 6-4 


RATIOS OF STANDARD DEVIATIONS OF HORIZONTAL VELOCITY 
COMPONENTS TO FRICTION VELOCITY; A FOR 
LONGITUDINAL, B FOR TRANSVERSE 


Time (min) A B 
1 Ocal5 O72 
3 0.20 0.17 
10 0.28 0.25 
30 0.39 0.37 


As expected, the open ocean values are much less than 
those over flat terrain, about one order of magnitude small- 
er. Since A and B deviate from the Panofsky and Dutton 
(1984) values, there is a strong possiblity that other pro- 
duction mechanisms independent of zZ, must exist. Possibili- 
ties could include buoyancy and mesoscale effects as 
proposed by Schacher et al. (1986), or terrain/orographic 


effects. 
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By assuming that the BLM data were obtained in the 


surface layer over uniform terrain, Eq. (24) becomes 


Oy/ue = Oy/Ux , (25) 


assuming A is about the same size as B. 
Using Equation 25, a final comparison was made using the 
ratio of o)/cy- The values used in Panofsky and Dutton 


(1984) resulted in a ratio of 
Oy/oy = 1.25 + 0.05 . 


The ratios 6,/oy are significant because they address the 
actual growth of a puff. 
Table 6-5 shows the q,/aq, ratios for each time average 


based on the processed BLM data. 


TABLE 6-5 


RATIOS OF STANDARD DEVIATIONS OF HORIZONTAL 
VELOCITY COMPONENTS 


Time (min) Oy/ Wy +0.05 
iL oe. 
3 Aero! 
10 1.24 
30 LS PAG 


As shown in Table 6-5, the o/o, ratios based on the BLM 
data are about the same as Panofsky and Dutton (1984) for 


overland. Note that the values in Table 6-5 decrease as the 


12 


time average increases. Based on Panofsky and Dutton 
(1984), it is expected that the 9,/y, ratios for the 30- 
minute time averaged values should be approximately 1.25. 
Using the parameterization constants from Table 6-2, %/%y 
ratios for the various production mechanisms can be 
computed. Shown below are the ratios of each production 
mechanism in the longitudinal (u) vs. transverse (v) 
directions (u/v) based on the empirically derived regression 


variables for the 30-minute time average. 


Shear = 105% 
Buoyancy = 1.41 , 
Mesoscale = 1.08 . 


The 06 ,/°%, ratios shown above indicate consistency with 
previous results which indicated that puff growth due to 
various production mechanisms are greater in the downwind 


direction versus the crosswind direction. 


C. RELATIONSHIP TO OTHER WORKS 

It was found during the course of fitting the data that 
the results are very similar to those found by Schacher et 
al. (1986). In all three cases (Ou v6), the short time 
average results are strongly dependent on the value of C, 
while the long time averages are not. Both wy ranges can 
not be fit without including shear, C,,. The conclusion is 
that shear production is not important for large averaging 
times since large scale processes dominate. Both the shear 


and buoyancy terms in Eq. (23) are functions of wind speed. 


74 


Shear production dominates at high wind speeds. As 
demonstrated previously in Table 6-1, the shear production 
contribution is greater than either buoyancy or mesoscale 
contributions at high wind speeds (>11 m/sec). 

When using the results of this thesis to parameterize 
Wind variablity one must first determine if the situation is 
stationary. In this context, stationarity means that there 
is a well Berecl ished wind, that one is not in a 
land-sea-breeze transition period. If stationarity exists 
Cns = 0, and if non-stationarity exists, the full Equation 
23 should be used to calculate the appropriate o. 

Due to the nature of the collection methods, sampling 
location, and local conditions, these findings may be less 
than precise. Therefore, certain caveats are appropriate: 


1. These results are only applicable to the coastal, 
overwater regime. 


2. These results may be location specific since the data 
were from two cruises but at a single distance 
offshore. 

3. The mesoscale production term was obtained by 
utilizing data obtained when the flow was driven by a 
sea-breeze cycle. This term may not be applicable to 
other conditions. 


4. The mesoscale term must be set to zero, Cypne = 0, for 
stationary conditions. 


Based on Schacher et al. (1986), there exists reasonable 
confidence that the buoyancy and shear production terms are 
transportable to other conditions and locations. There may 
be some difficulty with the buoyancy term since there exists 


the potential that not all the mesoscale influence was 


fos 


absent when the values were determined but, for its current 
potential uses, there should be no significant error. The 
mesoscale term was a large effect which may be site 
specific; however, that limitation does not preclude the 
possibility that this effect may extend far out to sea. 

The findings of this thesis are related to the results 
from Davidson (1974), where turbulence data were obtained 
over ocean waves during the BOMEX experiment of 1969, 
located 200 miles east-northeast of Barbados. Surface layer 
wind fluctuation and momentum transfer were shown to have 
been significantly influenced by both stability and wind 
wave coupling. The low frequency turbulence described by 
Davidson over the open ocean existed at low wind speeds and 
could have been due to any of several sources, including 
mesoscale production mechanisms. The similarity between the 
two sets of results is that low frequency turbulence at low 
wind speeds was observed in both near coast (Schacher et al. 
(1986)) and open ocean regimes (Davidson 1974). The 
implication is that this phenomenon may be generated by 
mesoscale production mechanisms, sensible heat flux, 
Wind-wave coupling, Kelvin waves, or a combination of 
sources. The recommendation then, is to conduct further 
research to establish continuity between the findings of 


Davidson (1974), Schacher et al. (1986), and this study. 
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D. PUFF MODELING IMPLICATIONS 

A final recommendation is to provide this work to the 
Naval Environmental Prediction Research Facility in 
Monterey, CA for incorporation into its Chemical Weapon 
Hazard Forecast (CHEMFO) model. Appendix F contains 
recommendations for the application of this thesis to modify 


CHEMFO such that it includes rudimentary puff modeling. 
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APPENDIX A 


SCATTERPLOTS 


The following are scatterplots of BIM data for 4, ,/U 


versus U. 
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APPENDIX B 


AVERAGED PLOTS, UNRESTRICTED CONDITIONS 


The following are averaged BIM-3 and 4 data for 


unrestricted conditions. 
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APPENDIX C 


AVERAGED PLOTS, STABLE CONDITIONS 


The following are averaged BLM data during only stable 


conditions. 
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APPENDIX D 


AVERAGED PLOTS, STATIONARY CONDITIONS 


The following are averaged BLM data plots during only 


stationary conditions. 


ry 
= . 


a 


1 min Reve 
Total 
STATIONARY 


ai 


Sigma-U 7 U 





WIND SPEED (m/sec) 


U % SIGMA Std Dev 

9 Q @.000 0.000 
ies Q @.000 @.200 
2e9 2 O41 882 
oo 7 O41 001 
ee ae 845 801 
Sieg 20 O41 @.000 
Bao od 038 @.800 
ted 38 835 @.00@ 
9.@ ib 041 @.80@ 
11.0 15 045 801 
15.9 ie 047 O01 


110 


3 min Ave 


Total 
STAT LONARY 

5 ec 
\ 
= 
a 
= 
B) 
uw] 

G6 





G oat 18 ge 
WIND SPEED (m/sec} 


U t SIGMA Std Dev 

5 Q 2.022 2.2200 
fies Q 2.220 2.222 
2.5 2 063. 204 
3.5 % .@52 201 
4.5 17 061 201 
S.5 20 .@S5 201 
6.5 39 .049 @.20@ 
7.5 30 044 @.200 
9.2 27 .@52 . 201 
11.2 15 054 201 
13.5 13 053 201 


ame 


Sigman a 


1@ min Ave 
Porta 
STATIONARY 





Y 3 1Q 
WIND SPEED (m/sec) 


U ' $¥ SIGMA Std Dev 

5 Q @.000 Q@.000 
1.5 Q @.000 0.802 
2.5 2 ml oZ .018 
3.5 7 069 802 
4.5 17 087 882 
5.5 20 = LF 001 
6.5 63 G62 .001 
an 3@ OS O01 
9.2 Za 058 . 901 
11.@ 1S .063 201 
5 


13 .Q56 . 201 


112 


15 


Sigma-U 7 U 


—£ ab 


O~-~ ao AI nom Sf OLIN — 
MWVeae nnnunuanwan 


. 


30 min Ave 
Total 
STATIONARY 





LZ 


“-AIND SPEED (m/sec) 


Sais 


Q. 
Q. 


MA 


000 
000 


Safes 
101 
ehlsS 
.089 
wots 
. 966 
67 
070 
862 


113 


Std Dev 


@.200 
QB. 800 
833 
80S 
004 
982 
901 
801 
001 
901 
801 


iS 


Sigma-V 7 U 


l1 min Ave 
Total 
STATIONARY 





Y 3 1g 
WIND SPEED (m/sec) 


U 7 z SIGMA Std Dev 

oo @ @.000 ®@.0@0 
5 Q @.002 ®@. 202 
2.0 ze » O51 203 
oo u O21 ®@.2002 
4.5 17 O27 @.000 
So 20 O24 ®@.000 
5.0 og 024 @. 202 
ZS 30 .@24 ®@.2002 
ie 27 824 @. 000 
11.@ 15 026 @.000 
Wee Le 2025 @.000 


114 


Cz 


Sigma-V 7 U 


———_- atid 


W- a NOM PON — 
Mma omwnannan mM 


ae Se aes 


S16 


Q. 
@. 


MA 


222 
222 


Oe! 
.@26 
842 
034 
iOS2 
33 
932 
Bi 2 7 
O31 


115 


3 min Ave 
Total 
STATIONARY 





1@ 


(m/sec) 


Std Dev 


S 


Sa0oaodc sd s 


. 200 
. 222 
885 
. 022 
. 022 
. 222 
. 020 
. 202 
. 222 
892 
. 020 


ies 


Sigma-V 7 U 


16 min Ave 
Total 
STATIONARY 





Q } 5 | 12 15 
WIND SPEED (m/sec3 


U € SIGMA Std Dev 

25 Q @.220 Q.222 
lo @ @.200 2.000 
2.2 2 - 180 O11 
S20 i 849 001 
4.5 ti »@62 201 
Sao 29 .@45 @. 200 
5.5 33 943 @. 802 
TAS 39 ~.@43 @.099 
2.0 27 048 @.000 
11.2 1S 938 Q.2090 
1635 io 935 8.822 


116 


Sitgmea-V ~ U 


— ah 


A- ON DOM PON — 
moaooununnununwanm Mm 


-“ 


fj — 
8euMghN OG GB 


— NN Ol 
ut ~1 ~) © 


12 


3H min Ave 
Total 
Sil Ginel Oita 





1g 


WIND SPEED (m/sec) 


Sit 


@. 
Q. 


MA 


200 
800 


ges 
O55 
0390 
og 
657 
O54 
One 
0oZ 
939 


ey 


Std Dev 


Q. 


g 


QO@ 
800 
O31 
OQ 1 
O02 
801 
201 
801 
201 
021 
880 


S 


APPENDIX E 


FITTING CURVES 


The following are empirical fitting curves for stable, 


stationary and unrestricted conditions for BLM data. 
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APPENDIX F 


PUFF MODELING IMPLICATIONS 


The Chemical Weapon Hazard Forecasting (CHEMFO) model 
developed by the Naval Environmental Predition Research 
Facility in Monterey, CA is a puff dispersion model used to 
predict the spread of chemical warfare agents, and is based 
on the work of Skupniewicz et al. (1984b). This model 
assumes a surface release with no vertical limit to the puff 
spread. However, that puff model only uses comparisons of 
relative dispersion data sets to field meteorological data 
for its parameterization. This thesis work can be extended 
to puff releases, which have operational significance to the 
U.S. Navy. 

Skupniewicz et al. (1984a) determined estimates of short 
range diffusion in the horizontal plane for the transverse 


direction only, using 


oy(T) = oy *T * £y(T/ty) , 


where: 


the standard deviation of the cross-wind 
velocity component, 


Ov 


rH 
ll 


the diffusion averaging time, 
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fy (T/tyz) = a universal function, and 


ty, = the Lagrangian time scale. 
The function fy(T/ty) can be related to the puff length 
scale by 


T= (L/U) , 


where 


L = the dimension of the puff, and 


U = the mean wind speed during the time average. 


Sheih (1981) experimentally determined the "universal" 
function for various Pasquill-Gifford stability categories 


using 
fr eee eG Zee) 


where 
t'; = an "apparent" integral time scale. 


It follows that for short range diffusion in the 


horizontal downwind direction 
Ox(T) = Oy * T * f,.(T/tyz) : 


where 


ee 


OQ, = the standard deviation of the downwind 
velocity component. 

By knowing the aowntnd characteristic of puff growth 
based on %,(T), the Gaussian plume dispersion model may be 
improved. Unfortunately, the above argument is theoretical 
and cannot be implemented bedause the Lagrangian time scale 
cannot be quantified; therefore, fy, and fy are not known. 
The alternative is to utilize the relationship between 9%, 
and O% as determined in this work and project’ the 
relationship between the needed 4% and °%, used in puff 
modeling. 

Skupniewicz (1984a) determined only transverse values of 
Oy based on tracer experiments where actual meteorological 
measurements were made during selected meteorological 
conditions. This thesis presents a regression equation to 
parameterize G@, and 0, under many different meteorological 
conditions. By assuming that the relationship of turbulence 
to plume growth is the same for both the x and y directions, 
the scaling of o, to Oy is expected to be the same for all 
Oy. This expectation. is also applied to the scaling of oy 
to Oo, for all 9. A relationship between o, and o, is 
presented earlier in this work and may provide an analogous 
relationship bewteen oy and ox. 

This thesis has dealt with puff growth using au, v grid 
based in a relative coordinate system where u is the axis 


along the mean wind flow. Based on this thesis, and 
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Panofsky and Dutton (1984), the ratio of 0,/ 5, is 
approximately 1.25 +0.05. Since Oy has been established as 
a function of oy, and y, is proposed to be a function of oy, 
it follows that there should exist a relationship such 
that oy/ oy is also approximately 1.25 +0.05. 

Incorporation of this relationship into the CHEMFO model 
should provide an initial attempt to more accurately portray 


puff growth over an open ocean. 


2 


LIST OF REFERENCES 


Businger, J.A., 1973: Turbulent Transfer in the Atmospheric 
Surface Layer. Workshop on Micrometeorology, AMS, 67. 


Davidson, K.L., 1974: Observational Results on the 
Influence of Stability and Wind-wave Coupling on 
Momentum Transfer and Turbulent Fluctuations over Ocean 
Waves. Boundary Layer Meteorology, 6, 305- 331. 


Garratt, J.R., 1977: Review of Drag Coeffiecients over 
Oceans and Continents. Monthly Weather Review, 105, 
915. 


Hojstrup, J., 1981: A Simple Model for the Adjustment of 
Velocity Spectra in Unstable Conditions Downstream of an 
Abrupt Change in Roughness and Heat Flux. Boundary 
Layer Meteorology, 21, 341-356. 


, 1982: Velocity Spectra in the Unstable Planetary 
Boundary Layer. J. Atmos. Sci., 39, 2239. 


Kaimal, J.C., J.C. Wyngaard, D.A. Haugen, O.R. Cote, Y. 
Izumi, S.J. Caughey and C.J. Readings, 1976: Turbulence 
Structure in the Convective Boundary Layer. J. Atmos. 
Sci., 33, 2152-2169. 


Kaimal, J.C., 1978: Horizontal Velocity Spectra in an 
Unstable Surface Layer. Atmos. Sci., 35, 18-23. 


Naval Postgraduate School, 1980: Naval Postgraduate School 
Shipboard and Aircraft Meteorological Equipment. 
Report 61-80-017 PR, Naval Postgraduate School, 
Monterey, CA, 22 pp. 


Panofsky, H.A. and J.A. Dutton, 1984: Atmospheric 
Turbulence, lst ed., John Wiley & Sons, Inc., 397 pp. 


Schacher, G.E., D.E. Spiel, cC.W. Fairall, K.L. Davidson, 
C.A. Leonard and C.H. Reheis, 1982: California Coastal 
Offshore Transport and Diffusion Experiments- 
Meteorological Conditions and Data. Report 61-82-007, 
Naval Postgraduate School, Monterey, CA, 380 pp. 


Schacher, G.E., 1986: Parameterization of Overwater 


Horizontal Wind Variability. Report 61-86-0114, Naval 
Postgraduate School, Monterey, CA, 211 pp. 


140 


Sheih, C., 1981: Pasquill-Taylor Dispersion Parameters 
Overwater Near Shore. Atmos. Environ., 15, 101-105. 


Skupniewicz, C.E. and G.E. Schacher, 1984: Measured Plume 


Dispersion Parameters Over Water, Volume I. Report 
61-84-012, Naval Postgraduate School, Monterey, CA, 101 
PP- 

, 1984: Assessment of the Performance of an In-Field 
Gaussian Plume/Puff Model for Overwater Use. Report 
61-85-002, Naval Postgraduate School, Monterey, CA, 55 
PP - 


Vandenberg AFB, DET 30, 2WS, 1982: Terminal Forecast 
Reference Notebook. Vandenberg AFB, CA, 98 pp. 


Aw 


INITIAL DISTRIBUTION LIST 
No. Copies 


Defense Technical Information Center 2 
Cameron Station 
Alexandria, Virginia 22304-6145 


Library, Code 0142 2 
Naval Postgraduate School 
Monterey, California 93943-5002 


Chairman (Code 68Mr) 1 
Department of Oceanography 

Naval Postgraduate School 

Monterey, California 93943-5004 


Chairman (Code 63Rd) a 
Department of Meteorology 

Naval Postgraduate School 

Monterey, California 93943-5004 


Director Naval Oceanography Division i: 
Naval Observatory 

34th and Massachusetts Avenue NW 

Washington, D.C. 20390 


Commander 1 
Naval Oceanography Command 

NSTL Station 

Bay St. Louis, Missouri 39522 


Commanding Officer 1 
Naval Oceanographic Office 

NSTL Station 

Bay St. Louis, Missouri 39522 


Commanding Officer J 


Fleet Numerical Oceanography Center 


Monterey, California 93940 


Commanding Officer 1 
Naval Ocean Research and Development 
Activity 


NSTL Station 
Bay St. Louis, Missouri 39522 


142 


10. 


Pl. 


12. 


i) 


14. 


15. 


16. 


17. 


18. 


19. 


Commanding Officer 

Naval Environmental Prediction 
Research Facility 

Monterey, California 93940 


Chairman, Oceanography Department 
U.S. Naval Academy 
Annapolis, Maryland 21402 


Chief of Naval Research 
800 N. Quincy Street 
Arlington, Virginia 22217 


Office of Naval Research (Code 420) 

Naval Ocean Research and Development 
Activity 

800 N. Quincy Street 

Arlington, Virginia 22217 


Scientific Liaison Office 

Office of Naval Research 

Scripps Institution of Oceanography 
La Jolla, California 92037 


Library 

Scripps Institution of Oceanography 
P.O. Box 2367 

La Jolla, California 92037 


Library 

Department of Oceanography 
University of Washington 
Seattle, Washington 98105 


Library 

CICESE 

P.O. Box 4803 

San Ysidro, California 92073 


Library 

School of Oceanography 
Oregon State University 
Corvallis, Oregon 97331 


Commander 

Oceanographic Systems Pacific 
Box 1390 

Pearl Harbor, Hawaii 96860 


143 


20. 


oak 


22 © 


23. 


24. 


Ze 


26. 


Zoe 


28. 


Zo. 


Chief, Ocean Services Division 

National Oceanic and Atmospheric 
Administration 

8060 Thirteenth Street 

Silver Springs, Maryland 20910 


Library 

Moss Landing Marine Lab 
California State Colleges 
Sandholdt Road 

Moss Landing, California 95039 


Library Acquisitions 

National Center for Atmospheric Research 
P.O. Box 3000 

Boulder, Colorado 80307 


NAVSEASYSCOM 

Damage Control Branch 

Code 55X25 

Bldg NC-2 

Washington, D.C. 20362-5101 


Professor G.E. Schacher 

Dean of Science 

Naval Postgraduate School 
Monterey, California 93943-5004 


Professor K.L. Davidson, Code 63Ds 
Department of Meteorology 

Naval Postgraduate School 
Monterey, California 93943-5004 


Mr. Tom Yencha, Code G51 
Naval Surface Weapons Center 
Dahlgren, Virginia 22448 


CDR James D. Branum 
2682 Olivestone Way 
San Jose, California 95132 


LT Tim Dowding 
6 Robin Way 
Westerly, Rhode Island 02891 


Dr. Donald L. Shearer 

TRC Environmental Consultants, Inc. 
8775 E. Orchard Road, Suite 816 
Englewood, Colorado 80111 


144 


20. 


31. 


S2. 


eS 6 


34. 


Dr. Bert Galloway, Code 3917 
Naval Weapons Center 
China Lake, California 93555 


Dr. Gloria Patton 

Theater Nuclear Warfare Project Office 
Department of Navy 

Washington, D.C. 20360 


CDR S.G. Colgan, Code 420B 
Office of Naval Research 
800 N. Quincy Street 
Arlington, Virginia 22217 


Mr. Dirk Herkhof 

Minerals Management Service 
Pacific ocsS Office 

1340 West 6th Street Rm. 200 
Los Angeles, California 90017 


Mr. Robert Harrison 

Western Oil and Gas Assn. 

727 West 7th Street 

Los Angeles, California 90017 


145 














2 
3 





°—! bs ates al ™~ @cr% re. 
br eamghiee i eee Oo ey ee eee 





4 i by 
Py-atealeeata ee DO e ee ee ts Se od 
y engunienmencamadel e eT 2 oh BALE 
Liege tienen ak 





tied ie pein | ae 
jure Dinetueedencel initia Le a i a 
hp-td ie Cate Ts Ae RPOADAS Mes 





thesD70153 





0 
Tras | ike mi F Parameterization of horizontal wind vel 
a TRI rer ecp eter ar eae ee wees ’ | ! 


ipinpketterdine LL et se) 


nihil | 
Oe er re t VMTN nn | 1 Wh \| | | 
rary pe ete at eT WW AMAT 

el waiit re roa ve oe Sank | | } \ | | iH | 
ees Perea ree 4s pee. EN | HI | | | 

a Sapte Sabet ry Tr 

Fel apeib ot 


whee oe eS 
EPPA LO ott ont 


2768 000 77213 1 
~ UDLEY KNOX LIBRARY 





Satapepadinn ae 
Meee te a 
ae Td See 


eihagheete aie eta ee A eT. 

a ee de eed ehelid ear mabtao, ee ns te 
ett teem Catt Te 6 « Ce a Mee | a Fir} 

on! aN en edt ee ey 


ed ae 
i ee Roh ome e686 del 1h | 


salt end en chon ee oe ad ~ 





1) rT A 





mI . 
ttn: Ce eee Py er 7 

pt gabe teats a en isete le ech ol 

oudeommee tO ee 





epee tLe Ye Eee 
anttinnahaise tment ee Saetrvedar hk ee ee ee 


L 
Siegen Semaine Lt ee ee ee | © chee h © Hoek ofan 

masta tenia nan hn tt ee eet ee eg ee ee ee hed.n "ge ee 
nk. Kemetic BP OI2 AD 991 ue Dee ee me ate, a oe a oe | 

cette eh dt ee ie hae kL het See TT ee eee Cy Pe | 


AP POE AD She on ee nl) oe med Or i fe te | 
7 cn 1a Pee Bn A 

















rd Lett Lit er 
eee ited 1 ed  hlleat at t . 
Ate anette tT cn oe ee er er 
ttn auteretdem ena ee ghana eet ee ee ae 
aioe Le 160% 0 of FEE. 0 oD epee 
iiiiedineenied ae sit-epgeaenel ee ee 
stile ae ep peree diate nai ahi Tt eee het be 
neu l-qteectenad ate oe om or: nme ee ah om.ee 
oobisimetinodans al tree cad) Lad $a! esometIncyee’ ea 4 
threaten 1s at Te ot ee ee oe Pot et at py eee 
wi teenaent oe tS LE teeth at ee 

iain ontest oe Sa dle ot aS 


ee er ae 











*.% hl ospda’s Ly _ 
5° Off, dan bile se? ae Lee ee 
ettepeied te eed Pt ht LEE eee ey? Bas # fe myt ng 
sdaelilate inten” Te ti 4 be eee to ee ee ee 7 od rc 
lr teblerstetndiinnd hin Le Pe y pry eran pe ee tert ae cae Ce ys 
en eth Let AY Lee dh en BAVA Mol e1@er, Lat XY Pal 
Syietllapedee dames cena) tet et. ee ee Ce tart bee LT) 4.1 Moy idet es 
eadopeginecoaahadet eee eee ee lie ie Let re 1a en wat os blah) 
epitaplilnm amin od ot Pee PP ee ob ed Le 1 Dw eo O48 hee te 

a belie ee et ie ee te a Pe ee ie ee oe Pe ed Cth ae 
petition apne be eppdina aedededer tal ae Le LE a ey 
oe) ee ete Lete® Jats Le ee kbd {11 | Sd Ot Oe Te) a] 
imine te ee beable be etnias aed, Ae LY Pere ooo peed Tao yt | 
ctebntittivemd Lh sett oe penta Mensa hal on oll Lt ee ee le te | are oe ee pm: 
pein ieiaiabgimae con en ee Ce eat ale een taeda et ae aot 
ee an Tre aS €26GBMA® 044, ot Oh wah wee i Eee ee ee ay 

hr Tea gilaheetheteee de ee ee Se ey es 
edie eee alee Eh ee 
shadinte papel eee ee 


i eae ey 
Slee eee ed 
























Par 
Led ea eo | 
homed dl od = PL Ee ae 
ee ey 
- Pa «Seale orm ey 
bath ieee Le Ae eo ee 
a (leet an tac ee, COR ngttinm ike okt Li, Ce ee 
eitietanale ee a hada dealetietden ne 1 OO eS Le ye ee Py 
sitirpmcdt nha dodo a statin I. a Yl | cP Ley eee ee ar ae | 
iubhpuntintee ait ie ae ehtel ait) tt ee ee ar 
peda 2 PRR Ree Oh Mm ones : Te « 


dor sherstcrel ate I Te Pe as _ at Pe eer De 
alien ee J 


- ae Tee et Te 
aed ee ae ed * #  MrBote mm BAG te mee 0 8 ood Py 
5 gh liaioan od cae rae a 1 ame pwd eas hg ra) 

bk dee Se tet ioe tee BaeGe Ss F yp want 
ail cinel atta Pibtimipdenblln ak dt eT a ee abled lamer tl oe aa 
mepuapdrpbhineladin i OY ot et eer “eile ice et oe eee ee 
at te eee eisai ot ee 


ar . wed gle eT ae tee 
edie oe fo teen eed om be iba 


ee er a rr ee a ak et ae oe 
Oe ee eS pee LY te ete ee ee San *s 
serine. ee ne 5 id 

Rate ob te  weatiatlahocell ott) 

nOtb that nh re wien 

— ee eee re  ) 
at te 
eed 
Sar he 











. dng a OOD 


































eee ee a ee | 










_— 





' 
ed 


er 

















ala 2 






i) Sie wnf otal Mh 
ee i tye 
et ee ee 
wh ptt g 








tet aie’ ta? aed 
eg ne eee 
ae ee 
~ = seirueiete ee 
Sea 





: et Peep 
ie ee ere err. 
ated ok ee ad 
sililn eat ead ac 
aed ek ea 

et ee ee 


steniehaliaat al fe a 
4 ’ er 
Ot 





| beet Bary ye 


sttetneell Anat Fant ot 
Et ta ee eas is 
ln ee we ee alae Ds Tl pitted oe 
id WA es ae > Toston he, Yee Bente iy he ies, 
beanie eT De eT 
> Soa ag, 7. phd te te ee Pn) bt ate 
hile ede Lo nt i a eS fy Mos ay 


"35. yey >. tt tht tah Bole Te ten LiL a ee *ofem ergs 
bh dtnash te te kitten ke ek eT ae 


aphnptealibehi inte ta 
MS ee Ryd bah teeta da hes artes Laat edes TU ae es 
sede ba de LT ree 

















Co 
* 907% e'hre we ow cs: 2 
Lie Le tary t oe wg% e6, 1 ¢ Wy 
oily seem & pl 
sw AVats nm gy 1¢ > 


raw 













i ahaetiecetmele to ee 


= bol priatao hate? aaltunlaht on sania a te A re 
b eee or We deka Lh. ton in SE 
bee ted hn tet ot 


= 5 . 5 , f 
AIS, OP “ ° miei n ‘ Py 
bed ta Ce et Cre oI y 4 e 
bet ar * 6 50 es -™ - 
fot th ee ee 
ey Ts 9 Pe dh tl Le Doc 2 #9 +0 A 
Peso ot) peer air Le a kD As fo rue @ F 
ofaye a’ 8 COLL Y | ® mt 
r he ie od at 
i ber. Slee ey Pt 
et Awe 
Fetreen ats 
ae 





Si Tn 
+ ~Pe teres Seg mr 
baad Tete | hed eee ee Let ey 
at Témy Pee P ae 
el ie tire Sa te te Th  Mpembe 4:5, » 
om oy! oft RK Rie es te a a 
tet nh eH he he 2 ae 9 | Sees vey 
bith il Moelle tit ee Te 8 int dhe Balin Teh Cs tS et Od -% 
plik te de te tl, Eee} Faded Soa te  , Y ee ae mi tne a: 
Painter td tepth te hd tet te batik i fed abet betadetn os A Se se? 
Se ee litle en “winraheriDe a eumr: RH SOme ome ity he art é¢ 
ipnditatiiem hs Lt Let he he th ie Soe Ltete ry 
eth tee th a tele hee fn ee itledches eo he ee TER RS 
mapdstintlihecadgunke tn en ee beta len els 


‘~~ OPE TR vice 
it od be Se hE 
E dh Lente oe et LL 
bocaptbdibatehendn bs mat et - 


reheat to Lae Ta 
oP a 04 


bed Be et tear 
aL iri bl ere 
PO owas 







os 















" 1 
a 2! e J 

Delt taal : L] way La 1 n 

ae re biged 0 0 tk 

el el ee 

Hele ee ia) 

a rh rl a arate 

he et ee 

. O@xse aletig © 

he Pf pin : A . ed 

4 <The a4 r f o 
eh a te 4 f Ca n D ‘ o ' 

ee) hes: ye it r a a} a” 

HRB. yhoo eg Sint: wae OFe i 

nad ih deed TE) k ST air ‘ e 

ala Lb aes FL er ee J ar . or i ry = @e ' 

date) ot oy borers ¢ A. ri nh Yet Ta | LI Dat og ' Lr ae | i 

in, Kl tet Abo tt Tee Sheers a aT) s0td “tw gt ' ws Cn er c 

Bits bret 2 a | a % 























St ey 
Swine tii cR woe 
































on Shh dae 


a 
hte desks eeden 


eee py 
bn etal 













i ee Thy 
Rete i hk Sh 
Lele ie kota to Kn 
sed h tek Lie LEE I 
@ eal! O¥etels ye a awry a ke Pa Sh 
Mediate i Midi Sap betel ed EY hm © pl ye, lupe l yee 8 
dt ae he Ths tT) | 
tre eer i 

Fons "CaS O26 poets ie | 
Vere a's won / "1 
A Rie Uae 1 8 wen 
aha Wrap wr Ty LP , 
bts Made tt ee a Rd 0 ey ea 
eee eS ed 












™ shew ive 
We uibvect mre 
pede! hls et tet te he | 



























Se Tn = q 

bya dedie Lessa dilee eet ta tL ae Tet 

bdtacltpnedbe ede Lats he Te ie) 

Spdep init bn 4 ot td hath -deletditcheeh oo TT TTT Sts Fs ae 

ht eaten ho te he. bath tht ede ttel Lio ne te 
owiiecincicie Ain ten en tel Re) 





te 

he Dk TL re Ue) i 
ih. ee ee LP | ¢ 
tthe We ie tT ee ee oer) e cn 7 
2. he ti to Serre ew ofeo~w Pt we m] 5a ¢ 48 . Ler 
bee ds Me th ee oe er i er ee eee | eget 
bbs Dek” Peal Rog ty tr a suey ¢ i 

3 °h & @agis @ See ws 4% q i a eT | Oe 1 gree 
3 Phi eth ih ie. Ue Ce ee ee 
7 at Me) a el 




















| Rk 68 @ 

bite te Ld We Me Woke TL 
Nes Scent Bx a 

S1Pwor eer as elwege feeete om 

ee 

rote 





























Ce ae 
poh El eh a en te ed 
eit sa oe tae y] 
Loire cite. 
web éersteve 


bh hed Inet dicho LIT ok baw Bk Me on 





an Loe 
Lh Fed 
: a ry 
eve + e LY Ae 

Pi *? er re) 1} Ce 
RA a a i ee ee i ie et ee ry ry 
Seah Te Pohnall tent Mh Then badd Eo A ee Do ae eo bas dtp ay ba ee er oe re A " 
tere we J ee Ce ee \ ee Sl ee I Sor ta | P 

Lithia dee te ee LTT. r e@ireae s ote 

RRL EY fore a a Md L 
i tt et et oT) 
Ke ie Ll 
r 














A 








O'Hare ow awe: 








blink Lh sche te 
OC ey Eh) ws: 














































PO RD LEP t' 6) ‘Yrwrg & ae eee 
+4 Pabedan te te eT eet he eye 
dnd Seog dod Rae ne te Mth datid ial d'm alin, 
1MOTP ah etre sd biel dei ahh TALLEY Bh 
aida 








' J 
rs 
i vee ¥ Py wy xis a4 Fieger & . 

i. Ai ee eo ee ee 7 4 

°F V%FS & oe. Biew.c eee Mah A ie RL | 

a LL mak oe Se A " 
i th i 
ror De teteeany 
OR OG 807s Hyer 
A Modded bl didi Pe Pe 
hth fa WUT ee @ Solr, be ee ie er ee 
S’aerns LORI el tee pie mites be 
inn 4h ds de Kee 0) Die Corty-@’icara e* 
hth eka os | @ ove aH A ded 





ahs bl #er 
i ee Loh iY a 

ee he Le es ? 18 
A i a he ae a a r ty 

a] ee ' 
er Ly ee a ee UY 
Te Venom ema ees y o 1 6H ‘a! 
git ae LY a) a i 

eam A Or +: 4 . 
Se te 
wu 















4 ea Xd Rad ot ee 

ep tta ae teh TS Tat tet pled Ala 

hte and oY 

ata ey Chl 
A} a 













Ad De | 7 i ’ rr | 
L 










































































